Relativistic shock waves in viscous gluon matter 
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We solve the relativistic Riemann problem in viscous gluon matter employing a microscopic parton 
cascade. We demonstrate the transition from ideal to viscous shock waves by varying the shear 
viscosity to entropy density ratio rj/s from zero to infinity. We show that axi rj/s ratio larger than 
0.2 prevents the development of well-defined shock waves on timescales typical for ultrarelativistic 
heavy-ion collisions. Comparisons with viscous hydrodynamic calculations confirm our findings. 
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In the 1970's, shock waves were theoretically predicted 
to occur in collisions of heavy nuclei This phe- 

nomenon has been experimentally investigated 0] and 
subsequently observed Recently, jet quenching Q 
has been discovered in heavy-ion collisions at Brookhaven 
National Laboratory's Relativistic Heavy-Ion Collider 
(RHIC). In this context, very exciting jet-associated par- 
ticle correlations [1] have been observed, which indicates 
the formation of shock waves in form of Mach cones Q in- 
duced by supersonic partons moving through the quark- 
gluon plasma (QGP). If true, it could give a direct access 
to the equation of state of the QGP, because the Mach 
cone angle is given by a = arccos(cs/z;jet), where Cg is 
the velocity of sound of the QGP. The velocity of sound 
is related to the equation of state via — dP/de, where 
P is the pressure and e the energy density. 

Shock waves can form and propagate only if matter be- 
haves like a fluid. The large measured elliptic flow coeffi- 
cient V2 indicates that the QGP created at the RHIC 
could even be a nearly perfect fluid. This is confirmed 
by recent calculations within viscous hydrodynamics Q 
and microscopic transport theory Q which estimate the 
shear viscosity to entropy density ratio rj/s to be less 
than 0.4 in order to not spoil the agreement with the V2 
data. However, it is an important question whether the 
r]/s value deduced from V2 data is sufficiently small to 
allow for the formation of shock waves. 

In this Letter we make an effort to answer this ques- 
tion by considering the relativistic Riemann problem in 
viscous gluon matter. Using the BAMPS microscopic 
transport model (BAMPS denotes the Boltzmann ap- 
proach of multiparton scatterings) pJ3| we demonstrate 
the transition from ideal shock waves with zero width, to 
viscous shock waves with nonzero width, to free diffusion 
by varying the shear viscosity to entropy density ratio 
ry/s from zero to infinity. We estimate the upper limit 
of the T]/s ratio, for which shocks can still be observed 
experimentally on the time scale of an ultrarelativistic 
heavy-ion collision. 

The initial condition for the relativistic Riemann prob- 
lem consists of two regions of thermodynamically equili- 



brated matter with different constant pressure separated 
by a membrane at z = 0, which is removed at t = 0. 
Matter is assumed to be homogeneous in the transverse 
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FIG. 1: (Color online) The solution of the Riemann problem. 
At t = 0, the pressure is Po ~ 5.43 GeVfm"^ for z < and 
P4 = 2.22 GeVfm"^ for z > 0. The upper panel shows the 
pressure and the lower panel the velocity at time f = 3.2 fm/c. 

(x, y) direction so that the further evolution of matter 
becomes (1 -|- l)-dimensional. For a perfect fluid, i.e., 
r] — 0, this problem has an analytical solution [llj, given 
by the solid curves in Fig. [TJ For a thermodynamically 
normal medium, there is a shock wave travelling into 
matter with the lower pressure (on the right) with ve- 
locity Ushock larger than the velocity of sound Cg- The 
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region of constant pressure behind the shock wave is the 
so-called shock plateau. Here, matter moves collectively 
with a constant velocity Wpiat as shown in the bottom 
panel of Fig. [TJ Simultaneously with the creation of the 
shock wave, a rarefaction wave travels with velocity Cs 
into matter with the larger pressure (on the left). For 
an ideal fluid, the solution of the Riemann problem is 
self-similar in the variable ( — z/t, i.e., the solution as a 
function of ( does not change with time. 

The velocity of the shock front is determined from the 
relativistic Rankine-Hugoniot-Taub equations fl^ and is 
given by [11 1 



"^shock 



{Ps - P4){e3 + P4) 
(64 + P3)(e3 - 64) 



1/2 



(1) 



where P3 and 63 (P4 and 64) are the pressure and energy 
density on the left (right) side of the shock front. For non- 
interacting ultrarelativistic gluon matter the equation of 
state is e = 3P. For the situation depicted in Fig. [2 
Ushock = 0.645. 

The extreme opposite of an ideal fluid (i.e., 77 = 
0) is a gas of free-streaming particles (i.e., rj = 00). 
The solution of the Riemann problem for free-streaming 
particles is given by the short-dashed lines in Fig. [U 
which can also be calculated analytically 13|. From 



— 7(1,0,0,1;), where v = 
2^-1/2^ The pressure is de- 



the solution for the single-particle distribution func- 
tion f(x,p), one computes the energy-momentum ten- 
sor T^"' = ^{(Pp/E)pf'p'' f{x,p), where p^' = {E,p) 
is the four-momentum. The energy density defined as 
e = Uf^T^^Ui, is the largest eigenvalue of T^. The eigen- 
vector is then the four velocity in Landau's defini- 
tion [14]. In our case u 
T°^/{e + T^^) and7= (l-w 

fined as F = -Af,^T^"'/3, where A^"' = g^"' - u^'u" and 
gfj.iy _ (1,-1,— 1,-1) is the metric tensor. For systems 
of massless particles T^^ is traceless and thus P = e/3. In 
the free-streaming case the characteristic structure of the 
solution of the Riemann problem for an ideal fluid is com- 
pletely washed out, and a clear distinction between the 
shock wave and the rarefaction fan is no longer possible. 

In the following, we study the influence of the rj/s ratio 
on the formation and evolution of shock waves by solving 
the Riemann problem with the parton cascade BAMPS 
[lo| . BAMPS is a microscopic transport model which 
solves the Boltzmann equation p'^ 9^/(2:, p) = C'{x,p) for 
on-shell gluons with the collision integral C{x,p). In this 
study, we consider only binary gluon scattering processes 
with an isotropic cross section. We remark that, although 
perturbative QCD (pQCD) favors small angles in binary 
gluon scatterings, the use of an isotropic cross section ef- 
fectively implements pQCD gluon bremsstrahlung, which 
has a broader angular distribution due to the suppression 
of soft coUinear radiation [l3| . 

In our calculations we use a constant rj/s value. In 
order to achieve this, we have to locally adjust the cross 



section a, since the particle density n is not constant. 
The shear viscosity 77 is given by 77 = 4e/(15i?''') [ist , 
where the transport collision rate R*"" — na*^" — 2na/3 
for isotropic scattering processes [l^ . n is calculated via 
n = Nf'Uf,, where TV = J{d^p/E)p'' f{x,p). We obtain 



V = -eAmfp : 





(2) 



where Amfp — ^/ {nu) is the gluon mean- free path. Binary 
collisions imply that we cannot maintain chemical equi- 
librium. In this case, in kinetic equilibrium the entropy 
density is calculated approximately via s = 4n — n In A, 
where A — n/rioq is the gluon fugacity measuring the 
deviation from chemical equilibrium. For a nonvanish- 
ing shear viscosity, the system will deviate from initial 
kinetic and chemical equilibrium during the evolution 
17| . Gluons are regarded as Boltzmann particles, so that 



the number density in thermal equilibrium is given by 
ncq = daT^ /i^^ with dc — 16 for gluons and T = e/(3n) 
denotes the temperature. Sending 77/s to zero the cross 
section will be unphysically large. However, it serves as 
a useful test of the parton cascade in the ideal hydrody- 
namical limit. 

In Fig. [Tl we show the results for various 77/s values 
as computed with BAMPS, demonstrating the gradual 
transition from the ideal hydrodynamical limit to free 
streaming of particles. Remarkably, the ideal solution is 
reproduced to very high precision by the BAMPS calcu- 
lation for a small 77/s = 0.001. In this sense, BAMPS 
can compare with state-of-the-art numerical algorithms 



used to solve the ideal hydrodynamical equations pj| . A 
larger 77/s value results in a finite transition layer where 
the quantities change smoothly rather than discontinu- 
ously as in the case of a perfect fluid. The width of the 
shock front is proportional to the shear viscosity [ll |. 

As seen in Fig. [U a nonzero viscosity smears the pres- 
sure and velocity profiles and impedes a clean separation 
of the shock front from the rarefaction fan. A criterion 
for a clear separation, and thus the observability of a 
shock wave, is the formation of a shock plateau, as in 
the ideal-fluid case. The formation of a shock plateau 
takes a certain amount of time, as demonstrated in Fig. 
[2] where we show the pressure and velocity profiles at 
different times for 77/s = 0.1. Formally, we define the 
time of formation of the shock plateau when the maxi- 
mum of the velocity distribution v{z) reaches the value 
v-piat of the ideal-fluid solution. From the bottom part of 
Fig. [21 we see that this happens at t = 3.2 fm/c. This 
agrees with what we obtained in the bottom part of Fig. 
[1] From this figure we also infer that, for 77/s > 0.1, 
a shock plateau has not yet developed &i t — 3.2 fm/c, 
whereas for rj/s < 0.1, it has already fully formed. 

In order to understand the timescale of formation of a 
shock wave, we define the quantity [l^ [l3, [l§| 
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mfp 
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where L = t (t'shock + Cs) is the width of the region 
bounded by the rarefaction wave traveUing to the left 
and the shock front moving to the right (in the ideal-fluid 
case). Note that Amfp is not constant; we approximate 
it by its maximum value which is assumed on the low- 
pressure side (the undisturbed medium) in front of the 
shock wave. The quantity K can be viewed as a "global" 
Knudsen number for the Riemann problem. K goes to 
zero at late times, which implies that the medium be- 
haves more and more like an ideal system. 
We find a scaling behavior for the solution of the Rie- 



1.00 
0.90 
0.80 
^ 0.70 

D. 

0.60 
0.50 
0.40 
0.20 

0.15 

^ 0.10 

0.05 

0.00 





initial condition 




t = 0.64fm/c 


\ 


t = 3.20 fm/c 


\ 


t 6 40 fm/c 


\ ; 




\ 




\ ; 




■ . \ 












'\ 




'■■ \ 


ri/s = 0.1 


': \ 




\ \ \ 












.'■ / 


\ \ 


/' ' ; 


\ 


/ 




/ I 


> 
1 

i ; 
i 


' ■■ 


i 

i '; 


/ I 


\ 







-6 -4 -2 2 4 6 
z(fm) 

FIG. 2: (Color online) The time evolution of the Riemann 
problem for 77/s = 0.1 and the initial condition of Fig. [T] 

mann problem: The pressure profile P(2;, 77/s)/Po (and 
also the velocity profile) is only a function of C = z/t 
and K, i.e., P{z,t;ri/s)/Po = F{(;K), for a given ra- 
tio P4/P0. The value of K at which the shock wave 
forms, called Kf, is universal for a given P^/Pq. With 
Eq. ([2]) we obtain Amfp = 10/ {3X4)?]/ s in the undisturbed 
medium with the lower pressure. Thus, using the obser- 
vation from Fig. [2] that the shock wave forms at i = 3.2 
fm/c for r]/s ~ 0.1, we find that Kf — 0.053. 

By inserting A,„fp = 10/(3r)77/s into Eq. the for- 
mation time of shock waves is given by 



10 



V 



3 Kf {vshock + Cs)T s 



(4) 



FigureElshows the relation (g]) with T = 350 MeV and for 
various initial pressure discontinuities P^/Pq. The differ- 
ence in slopes refiects the dependence oi Kf and Wshock on 
the ratio P^/Pq. For -q/ s = 0.2 a shock will not be visible 
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FIG. 3: (Color online) Formation time of shock waves as a 
function of rj/s. 



until 6—7.2 fm/c, which most likely exceeds the lifetime 
of the QGP at the RHIC. Vice versa, if shock phenomena 
are discovered at RHIC, this could be an indication that 
the QGP has a small r]/s ratio, probably smaller than 
0.2. For a more viscous QGP no shock waves and thus 
no Mach cones will be formed. In a relativistic heavy-ion 
collision, the temperature is decreasing during the expan- 
sion. Thus, according to Eq. ([4]), tf is even larger and 
shock waves may be even harder to observe. 

To confirm our results calculated with BAMPS, we 
solve the Riemann problem within relativistic dissipative 
fiuid dynamics. To this end, we use the Israel-Stewart 
(IS) equations [l^l (see Ref. [2lj for different approaches). 
We neglect heat conductivity and bulk viscosity (thus as- 
suming local chemical equilibrium). In 1+1 dimensions, 
the IS equations reduce to 



5tr°" -t- d,{vT°°)=-d,{vP + im), 

dtT°' + d,.{yT"'')=-d,{p + n), 



-ydfTT + ^vdzTi^—{-KNS - 7r) - 77 

2 



61 + D In 



(5) 
(6) 

,(7) 



where = d^u^ and D = u^^d^. The laboratory frame 
energy and momentum density are given by T"*' = (e + 
P + - (P ^), and = (e + P + ^)-i'^v, while the 
system is closed by a massless gluon equation of state. 
For vanishing shear viscous pressure, tt, Eqs. (O and ^ 
reduce to the equations of ideal hydrodynamics. 

The first-order theory of relativistic dissipative fluid 
dynamics deflnes the viscous pressure algebraically bv the 
relativistic Navier-Stokes value ttns = —(4/3)77^? (1 j. [2^ . 
In the IS equations, the local shear pressure relaxes to 
the Navier-Stokes value on a timescale comparable to the 
mean-free time between collisions, given by the relaxation 
time = 277/32, where /Ja = 3/(4P) 

The IS equations ([5])- ([7]) are solved numerically us- 
ing the viscous sharp and smooth transport algorithm 
(vSHASTA) approach Figure [4] shows comparisons 
between the results from BAMPS and vSHASTA calcu- 
lations for r]/ s = 0.01 and 0.2. We see perfect agreement 
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FIG. 4: (Color online) Same as Fig. [T] Results obtained from 
vSHASTA and BAMPS calculations are compared. 



for rj/s = 0.01, whereas for the larger value rj/s = 0.2 
small deviations in the region around the shock front 
and the rarefaction fan are found. The reason can be un- 
derstood considering the Knudsen number Kg = XmipO. 
The IS equations contain terms of second order in Kg 
24| . However, at the shock front, the macroscopic scale 
over which the fluid velocity varies is comparable to the 
microscopic scale Amfp, such that Kg ^ 1, and one would 
need higher powers of Kg in the IS equations to improve 
the description. Wherever Kg is large, the applicability 
of the IS theory is questionable. Microscopic transport 
calculations do not suffer from this drawback. 

In summary, employing the parton cascade BAMPS 
we have solved the relativistic Riemann problem. The 
transition from ideal-fluid behavior to free streaming is 
demonstrated. Numerical results from BAMPS agree 
well with those obtained from viscous hydrodynamical 
calculations. We found that the formation of shock waves 
in gluon matter with rj/s > 0.2 probably takes longer 
than the lifetime of the QGP at the RHIC. Whether 
Mach cones from jets can be observed within nu- 
clear collisions will be studied in the future within the 
BAMPS approach. 
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